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ABSTRACT 

PKS 0521—36 is an Active Galactic Nucleus (AGN) with uncertain classification. We investi¬ 
gate the properties of this source from radio to 7 rays. The broad emission lines in the optical 
and UV bands and steep radio spectrum indicate a possible classification as an intermediate 
object between broad-line radio galaxies (BLRG) and steep spectrum radio quasars (SSRQ). 
On pc-scales PKS 0521—36 shows a knotty structure similar to misaligned AGN. The core 
dominance and the 7-ray properties are similar to those estimated for other SSRQ and BLRG 
detected in 7 rays, suggesting an intermediate viewing angle with respect to the observer. In 
this context the flaring activity detected from this source by Fermi-LAT between 2010 June 
and 2012 February is very intriguing. We discuss the 7-ray emission of this source in the 
framework of the structured jet scenario, comparing the spectral energy distribution (SED) of 
the flaring state in 2010 June with that of a low state. We present three alternative models cor¬ 
responding to three different choices of the viewing angles 6 ^ = 6 °, 15°, and 20°. We obtain a 
good fit for the the first two cases, but the SED obtained with = 15° if observed at a small 
angle does not resemble that of a typical blazar since the synchrotron emission should domi¬ 
nate by a large factor 100) the inverse Compton component. This suggests that a viewing 
angle between 6° and 15° is preferred, with the rapid variability observed during 7-ray flares 
favouring a smaller angle. However, we cannot rule out that PKS 0521—36 is the misaligned 
counterpart of a synchrotron-dominated blazar. 

Key words: galaxies: active - galaxies: nuclei - galaxies: individual (PKS 0521-36) - galax¬ 
ies: quasars: general - gamma-rays: galaxies - gamma-rays: general 


1 INTRODUCTION 

Active galactic nuclei (AGN) are compact regions at the centre of a 
few percent of galaxies with a non-stellar emission overwhelming 
the thermal contribution of the entire galaxy. AGN include a super- 
massive black hole (SMBH) as central engine, whose strong gravi¬ 
tational potential pulls the surrounding materials inwards, forming 
a disc of hot plasma. In addition, gas clouds move in the potential 
well of the SMBH, producing optical and UV emission lines. The 


central region is surrounded by absorbing material in a flattened 
configuration, idealized as a toroidal shape, located at ~1-10 pc. 
In radio-loud objects there is the additional presence of relativis¬ 
tic jets, roughly perpendicular to the disc. According to the Unified 
Model dUrrv & Padovaniiri995h . AGN types are classified on the 
basis of the orientation of their jet with respect to the observer. 
Blazars, usually divided into flat spectrum radio quasars (FSRQ) 
and BL Lac objects, represent the fraction of AGN with their jet 
oriented at very small viewing ang le, which causes relativist ic aber¬ 
ration and emission amplification telandford & ReesI 19781) . Radio 
galaxies should be the misaligned parent population of blazars. On 
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the basis of their radio morphology and power, radio galaxies are 
classified as Fanaroff Riley type I (FR I) or Fanaroff Riley type 
II (FR II). Decelerating jets and kpc scale edge-darkened lobes 
are found in the weaker FR I radio galaxies, while relativistic jets 
and edge-bri ghtened radio lobes are found in the stronger FR II ra¬ 
dio galaxies ( IFanaroff & Rilevll974ll . According to the Unification 
scenario for radio-loud AGN, FR I radio galaxies correspond to BL 
Lac objects, and FR II radio galaxies are the parent population of 
FSRQ. 


_ ^PKS 0521—36 was first classified as an N galaxy l Bohone^D 

1965h. and then a BL Lac object jPanziger et al.l 1 19791 : 


Burbidge et al.llT98^ . However, this source shows br oad and vari¬ 
able nuclear emission lin es in optical and UV bands iUlrich et al.l 
I1981L IScarna et'^ll995ll with equivalent width (EW) much larger 
than 5 A (rest-frame), the threshold historically proposed to distin¬ 
guish between BL Lac objects and FSRQ (e.g.. lStickel et al.iri99ll : 
IStocke et'^ll99lh . This suggests that PKS 0521—36 is a misclas- 
sified AGN. For this source, a large-scale optical jet well aligned 
with the kiloparsec radio jet was observed, wit h a clear correspon- 
dence between the radio and optical structures dScarna et alJll99lj . 
and references therein). In radio and optical bands the jet resem - 
bles that of the nearby radio galaxy M 87 (e.g.. lSparks et al.ll 19941) . 
Ch andra detected the X-ra y counterpart of the innermost 2-arcsec 
jet teirkinshaw et alj 2002 h . 

In the 7 -ray energy range PKS 0521—36 was tentatively 
detected by the Energetic Gamma Ray Experiment Telescope 
(EGRET) on board the Compton Gamma-Ray Observatory in 
Phase 1 dLin et al.l Il995h. but it was no t included in the Third 
EGRET catalogue ( Hartman et al.l Il999h . No beaming effect is 
needed for the core brightness temperature estimated in ra- 
dio, consistent with the non-detection of superluminal motion 
dTingav & Edwardil2002h . By modeling the spectral energy dis- 
tribution ( S EP) a nd taking into account its radio characteristics, 
IPian et aP d 19961) derived a viewing angle of 30° with a bulk 
Lorentz factor of 1.2. Very Long Baseline Array (VLBA) image 
showed that the same Position Angle (PA) found on the parsec- 
scale jet is maintained, without an y significant bending , over 3 or¬ 
ders of magnitude of length scale dGiroletti et al.ll2004l) . This may 
be a further indication of a relativ ely large viewing ang le. Eor a 
plausible Lorentz factor of F = 5, iGiroletti et al.l d2004 ) derive a 
viewing angle 6 in the range 21 ° — 21°. 

On 2010 June 17, a strong 7 -ray flare from PKS 0521—36 was 
detected by the Large Area Telescope (LAT) on board t he Fermi 
Gamma-ray Space Telescope satellite diafrate et al.ll201oh . trigger¬ 
ing a Swift follow-up observation that c onfirmed the high activit y 
of the source in optical, UV, and X-rays dP’ Ammando et al .l 2010 h . 
In this paper we investigate the nature of this object and its emis¬ 
sion mechanisms by the analysis of multifrequency data collected 
from radio to 7 rays, focusing in particular on the 2010 June flaring 
activity. 

The paper is organized as follows: in Section 2, we report the 
LAT data analysis and results. In Section 3, we report the results of 
the new Swift and archival XMM-Newton and Chandra data analy¬ 
sis. Radio data collected by the VLBA, Very Large Array (VLA), 
University of Michigan Radio Astronomy Observatory (UMRAO), 
and Submillimeter Array (SMA) are presented in Section 4. In Sec¬ 
tion 5, we discuss the source properties from radio to 7 rays. In 
Section 6 , we discuss the modeling of the overall SEP and draw 
our conclusions in Section 7. 

Throughout the paper the quoted uncertainties are given at the 
la level, unless otherwise stated. The photon indices are param¬ 
eterized as dN/dE oc where Tj^ is the photon index at 
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Figure 1. LAT light curve of PKS 0521—36 for E > 100 MeV using a PL 
model with photon index fixed to T-,. = 2.45 in monthly time bins from 2008 
August 4 to 2012 August 4 (MJD 54682-56143). The downward arrow 
represents a 2cr upper limit. The dashed line indicates the average flux over 
the whole period. 


the different energy bands. We used a A cold dark matter cosmol- 
ogy with Hq = 71 km s“^ Mpc“^, = 0.73, and flrn = 0.27 
jKomatsu et alJlrOO^ . The corresponding luminosity distance at 2 
= 0.056 is dr = 246.8 Mpc and 1 arcsec corresponds to a projected 
size of 1.073 kpc. 


2 FERMI-L AT DATA 

The LAT on board the Fermi satellite is a 7 -ray telescope operat¬ 
ing from 20 MeV to > 300 GeV, with a large peak effective area 
(~ 8000 cm^ for 1 GeV photons), an energy resolution typically 
~10 per cent, and a field of view of about 2.4 sr with single-photon 
angular resolution (68 per cent containment radius) of 0 ? 6 at £ = 
1 GeV on-axis. Details about the LAT are given by I Atwood et al.l 

1 I2OO9I) . 

The LAT data reported in this paper were collected over 
the first 4 years of Fermi operation, from 2008 August 4 (MJD 
54682) to 2012 August 4 (MJD 56143) in the 0.1-100 GeV en¬ 
ergy range. During this time, the Fermi observatory operated al¬ 
most entirely in survey mode. The analysis was performed with 
the ScienceTools software package version v9r32p5 Q. Only 
events belonging to the ‘Source’ class were used. The time in¬ 
tervals when the rocking angle of the LAT was greater than 
52° were rejected. In addition, a cut on the zenith angle (< 
100°) was applied to reduce contamination from the Earth limb 
7 rays, which are produced by cosmic rays interacting with 
the upper atmosphere. The spectral analysis was performed with 
the instrument response functions P7REP_SOURCE_V15 using 
an unbinned maximum-likelihood method implemented in the 
tool gtlike. Isotropic (‘iso_source_v05.txt’) and Galactic diffuse 
emission (‘glLiem_v05_revLfit’) components were used to model 
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Figure 2. LAT light curve of PKS 0521—36 for E > 100 MeV using a PL 
model with photon index hxed to = 2.16 in 3-day (black circles) and 
1-day (red squares) time bins from 2010 June 4 to July 4 (MJD 55351- 
55381). 


Figure 3. The 0.1-100 GeV flux in units of 10“® ph cm“^ s"^ (top panel) 
and the photon index from a PL model (bottom panel) for PKS 0521—36 us¬ 
ing 1-month time-bins. The dashed line in both panels represents the mean 
value. 


the backgrountl^ The normalizations of both components were al¬ 
lowed to vary freely during the spectral fitting. 

We analysed a region of interest of 10° radius centred at the 
location of PKS 0521—36. We evaluated the significance of the 7 - 
ray signal from the source by means of the maximum-likelihood 
test statistic TS = 2x(logI/i - logLo), where L is the likelihood of 
the data given the model with (Li) or w ithout (Lp) a point s ource 
at the position of PKS 0521—36 (e.g.. iMattox et al.lll996h . The 
source model used in gtlike includes all o f the point sources 
from the third Fermi-hPH catalogue (3FGL; lAcero et"^ l2015l) 
that fall within 15° of the source. The spectra of these sources 
were parametrized by power-law (PL), log-parabola (LP), or ex¬ 
ponentially cut-off power-law model, as in the 3FGL catalogue. A 
first maximum-likelihood analysis was performed to remove from 
the model faint sources with fluxes lower than lxl 0 ~® ph cm~^ 
s“^ for avoiding possible problems of fit convergence. A sec¬ 
ond maximum-likelihood analysis was performed on the updated 
source model. In the fitting procedure, the normalization factors 
and the spectral parameters of the sources lying within 10° of PKS 
0521—36 were left as free parameters. For the sources located be¬ 
tween 10° and 15° from our target, we kept the normalization and 
the spectral parameters fixed to the values from the 3FGL cata¬ 
logue. 

The fit with a PL model, dN/dE oc {E/, to the data 
integrated over 48 months of Fermi operation (2008 August-2012 
August) in the 0.1-100 GeV energy range results in TS = 4945, 
with an integrated average flux of (13.9 ± 0.4)xl0“® ph cm~^ 
s“^, and a photon index F-, = 2.45 ± 0 . 02 . Figure [T] shows the 7 - 
ray light curve of the first 4 years of Fermi observations of PKS 
0521—36 built using 1-month time bins. For each time bin, the 
photon index was frozen to the value resulting from the likelihood 
analysis over the whole period. The systematic uncertainty in the 
flux is dominated by the s ystematic uncertainty in the effective area 
jAckermann et al.l 2012 d) . which amounts to 10 per cent below 100 

^ http://fermi.gsfc.nasa.gov/ssc/data/access/lat/Background 
Models.html 
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Figure 4. LAT light curve of PKS 0521—36 for E > 100 MeV using a PL 
model with photon index fixed to F-y = 2.33 in 3-day (black circles) and 
1-day (red squares) time bins from 2011 April 6 to May 5 (MJD 55657- 
55686). The downward arrow represents a 2a upper limit. 

MeV, decreasing linearly with the logarithm of energy to 5 per cent 
between 316 MeV and 10 GeV, and increasiM linearly with the 
logarithm of energy up to 15 per cent at 1 TeVQ 

PKS 0521—36 has been quite active since 2009 September. 
An increase of the 7 -ray flux was observed in 2010 January, and 
subsequently a significant flaring activity occurred in 2010 June. 
A second light curve focusing on the period of the highest activ¬ 
ity (2010 June 4-July 4; MJD 55351-55381) was built with 3- 
day and 1-day time bins (Fig. 1^. Leaving the photon index free 
to vary during 2010 June 4—July 4, the global fit results in TS = 901 
and a photon index r.y = 2.16 ± 0.05, indicating a harder-when- 

® http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html 
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Table 1. Unbinned likelihood spectral fit results. 


Time Period (MJD) 

Time Period (UT) 

PL 

r 

TSpl 

a 

LP 

ft 

TSlp 

54682-56143 

2008-08-04/2012-08-04 

2.45±0.02 

4945 

2.40±0.05 

0.07±0.02 

4962 

55351-55381 

2010-06-04/2010-07-04 

2.16±0.05 

901 

2.14±0.07 

0.02±0.02 

902 


brighter behaviour, similar to wh at was observed in FSRQs and 
low-synchrotron-peaked BL Lacs jAbdo et alj2010bll . The peak of 
the emission was observed on 2010 June 30 (MJD 55377), with 
a flux of (148 ± 18)xl0"® ph cm~^ s'^ in the 0.1-100 GeV 
energy range, a factor of ~11 higher than the average 7 -ray flux 
during 2008 August-2012 August. The corresponding apparent 
isotropic 7 -ray peak luminosity is 8.5x10"*® erg s“*. During the 
period 2010 June 4-July 4, the fit with an LP model, dN/dE oc 
[E/ iog(B/Bo)^ 0.1-100 GeV energy range results 
in TS = 902, with a spectral slope q = 2.14 ± 0.07 at the decorrela¬ 
tion energy Eq = 384 MeV, and a curvature parameter /3 = 0.02 ± 
0.02 (Table[TJ. We used a likelihood ratio test (LRT) to check a PL 
model (null hy pothesis) against an LP model (alternative hypothe¬ 
sis). Following iNolan et al.l j2012h . these values may be compared 
by defining the curvature test statistic TScurve=(TSLP - TSpl)- The 
LRT results in TScurve = L corresponding to a ~l(j difference, in¬ 
dicating no significant curvature in the 7 -ray spectrum during the 
flaring period. Considering the entire period 2008 August 4 - 2012 
August 4 the fit with an LP results in TS = 4962, with a spectral 
slope a = 2.40 ± 0.05 at the decorrelation energy Eo = 384 MeV, a 
curvature parameter j3 = 0.07 ± 0.02, and an average flux of (13.2 
± 0.8) X 10“® ph c m~^ s~* (TablefB, consistent with the results re¬ 
ported in the 3FGL jAcero et alj2^5h . The LRT results in TScurve 
= 17, corresponding to a ~4o- difference (for 1 degree of freedom). 
This indicates a significant curvature in the 7 -ray spectrum over the 
entire period, contrary to what is observed in shorter periods. 

In addition to the 3-day and 1-day light curves shown in Fig.|2 we 
also computed a light curve in 12 -hour bins during the period of 
brightest flux (2010 June 26-July 1). In this light curve (not shown) 
we note a significant increase from (49 ± 15)x 10“® ph cm~^ s~* 
to (193 ± 29)xl0“® ph cm“^ s“* between the first and second 
12-hr bin on 2010 June 29. The peak flux at 12-hr time-scale corre¬ 
sponds to an apparent isotropic 7 -ray luminosity of 1.1 x 10 "*® erg 

S-*. 

We investigated whether spectral changes are present during 
the period 2008 August-2012 August, using a PL model. We plot 
the photon index against the 7 -ray flux above 100 MeV estimated 
on a monthly time-scale. Unlikely as in Fig.[T] in which the photon 
indices are fixed, for each time bin the spectral parameters for PKS 
0521—36 and for all the sources within 10° from it were left free 
to vary. No obvious relation between flux and photon index was 
observed (Fig.O. 

Other high-activity periods of PKS 0521—36 have been ob¬ 
served in 2011 February and April, and 2012 January. In particular, 
in 2011 April the source was not detected for the first part of the 
month, and then a rapid increase of flux from (25 ± 13)xl0“® 
ph cm“^ s“* to (185 ± 29)xl0“® ph cm“^ s“* was observed 
between 2011 April 23 and 24 (Fig.0 with F.^ = 2.33, the value 
obtained leaving the photon index free to vary on monthly time- 
scales (see Fig.[3] bottom panel). As for Fig.[2 we investigated the 
12-hr light curve during the period of brightest flux (2011 April 
22-26). We note that an increase from (129 ± 32)x 10“® ph cm~^ 
s“* to (275 ± 56)x 10“® ph cm~^ s“* was observed between the 


first and second 12-hr bin on 2011 April 24. The peak flux at 12-hr 
time-scale corresponds to an apparent isotropic 7 -ray luminosity of 
L2x 10^® erg s“*. It is interesting to note that during the 2010 June 
and 2011 April flares the peaks show a moderately asymmetric pro¬ 
file (i.e. different rising and decaying times) on a daily time scale. 
In particular, the rapid increase observed in 2011 April suggests a 
cooling time longer than the light crossing time R/c related to fast 
injection of accelerated particles and a slower radiative cooling. 

Analyzing the LAT events with E > IQ GeV collected over 
2008 August-2012 August the fit results in TS = 89, a photon index 
r.^ = 2.71 ± 0.27, and a flux of (1.43 ± 0.36) x 10"*° ph cm“^ s"*. 
By means of the gtsrcprob tool, we estimated that the highest 
energy photon emitted from PKS 0521—36 (with probability > 90 
per cent of being associated with the source) was observed on 2010 
August 17 (MJD 55425), at a distance of 0?03 from the source and 
with an energy of 73 GeV. 


3 X-RAY, UV, AND OPTICAL DATA 


3.1 Swift observations 

The Swift satellite jGehrels et al.ll200^ performed nineteen obser¬ 
vations of PKS 0521—36 between 2005 May 26 and 2011 Septem¬ 
ber 5. In particular five observations were triggered by the 7 - 
ray activity observed in 2010 June. The observations were per- 
forme d with all three on-b oard instruments: the X-ray Telescope 
(XRT; [Burrows et ai]|2005L 0 2-lC) .0 keV), the UV/Optical Tele¬ 
scope (UVOT; Roming et al. 20051. 17 0-600 nm), and the Burst 
Alert Telescope (BAT; Barthelmv et al.ll200^ . 15-150 keV). 


3.1.1 Swift-BAT 

The hard X-ray flux of this source is below the sensitivity of 
the BAT instrument for the short exposure time of the individ¬ 
ual Swift observations. However, PKS 0521—36 is detected in the 
BAT 70-month catalogu^E generated from the all-sky survey in 
the time period 2004 November-2010 August. The data reduction 
and extraction procedure of the 8 -channel spectrum is described in 
[Baumgartner et all ( l2013h . The 14—195 keV spectrum is described 
by a power law with photon index of 1.83 ± 0.15 (x?ed = 2.87, 
5 d.o.f). The resulting 14—195 keV flux is (3.5 ± 0.3)xl0“** erg 


3.1.2 Swift-XRT 

The XRT data were processed with standard procedures 
(xrtpipeline vO . 12.6), filtering, and screening criteria by 
using the HEAsoft package (v6.13). The data were collected in 
photon counting mode in all the observations, and only XRT event 
grades 0-12 were selected. Source events were extracted from a 
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Table 2. Log and fitting results of Swift-XKT observations of PKS 0521—36. A power-law model with A^h fixed to the Galactic column density was used. 


Date 

(MJD) 

Date 

(UT) 

Net Exposure Time 
(s) 

Tx 

Flux (2-10 keV) 
(10~^^ erg cms“^) 

X?ed (d o-f-) 

53516 

2005-26-05 

900 

1.72±0.14 

1.34±0.12 

1.034(17) 

54498 

2008-02-02 

4832 

1.61 ±0.06 

1.32±0.06 

1.062 (88) 

54503/04 

2008-02-07/08 

4947 

1.59±0.06 

1.18±0.06 

0.918(64) 

54509 

2008-02-13 

2901 

1.55±0.09 

1.27±0.08 

0.869(41) 

55260 

2010-03-05 

1681 

1.67±0.12 

0.88±0.11 

1.082 (20) 

55263 

2010-03-08 

2405 

1.62±0.10 

1.20±0.10 

0.887 (35) 

55366 

2010-06-19 

2759 

1.62±0.08 

1.86±0.10 

1.018(42) 

55370 

2010-06-23 

3951 

1.57±0.07 

2.38±0.14 

1.020(57) 

55382 

2010-07-05 

2774 

1.64±0.09 

1.97±0.13 

1.089 (42) 

55386 

2010-07-09 

2904 

1.66±0.09 

2.06±0.15 

1.073 (37) 

55390 

2010-07-13 

2822 

1.60±0.08 

1.90±0.11 

0.992 (40) 

55626 

2011-03-06 

1126 

1.78±0.16 

1.30±0.18 

0.776 (13) 

55627 

2011-03-07 

2015 

1.75±0.10 

1.49±0.09 

0.899 (32) 

55661/02 

2011-04-09/10 

1868 

1.62±0.08 

1.34±0.09 

0.968 (34) 

55804/05 

2011-08-31/09-01 

1628 

1.79±0.12 

0.80±0.08 

1.117(17) 

55809 

2011-09-05 

1393 

1.67±0.12 

1.77±0.15 

0.928 (21) 


circular region with a radius of 20 pixels (1 pixel = 2.36 arcsec), 
while background events were extracted from a circular region with 
radius of 50 pixels away from the source region. Some observations 
showed a source count rate > 0.5 counts s~^; thus pile-up correc¬ 
tion was required. For those observations we extracted the source 
events from an annular region with an inner radius of 3 pixels (esti¬ 
mated by means of the PSF fitting technique) and an outer radius of 
30 pixels, while background events were extracted from an annular 
region centered on the source with radii of 70 and 120 pixels. An¬ 
cillary response files were generated with xrtmkar f , and account 
for different extraction regions, vignetting and PSF corrections. We 
used the spectral redistribution matrices in the calibration database 
maintained by FlEAS ARCfl The data collected during 2008 Febru¬ 
ary 7 and 8 , 2011 April 9 and 10, and 2011 August 31 and Septem¬ 
ber 1 were summed in order to have enough statistics to obtain a 
good spectral fit. 

We fit the spectra in the 0.3-10 keV energy range with an ab- 
sorbed power law u sing the photoelectric absorption model tbabs 
jWilms et al.ll 2000 l) . with a neutral hyd rogen column densit y fixed 
to its Galactic value (3.58x 10^° cm~^: [Liberia et al.l2005ll . Con¬ 
sidering the lower statistics with respect to the XMM-Newton and 
Chandra observations, more detailed spectral modelling was not 
performed with the XRT observations. The fit results are reported 
in Table|2 Swift-XKT observed the source during 2005-2011 with 
a 2-10 keV flux in the range (0.8-2.4)x 10“^^ erg cm“^ s“^. No 
significant change of the photon index was observed. During 2010 
mid-June, in the period of highest 7 -ray activity, an increase of the 
X-ray flux was observed, but not accompanied by significant hard¬ 
ening of the spectrum. Furthermore, by comparing the 2-10 keV 
flux observed by Swifl-XKF with those observed by Chandra (see 
Sect. l3.2| for details) we noted a significant increase of the flux (a 
factor of ~10) between 1999 and 2011 (Fig.|5ll. 

3.1.3 Swift-UVOT 

UVOT observed PKS 0521—36 in all its optical (v, b, and u) 
and UV (uvwl, uvm2, and uvw2) photometric filters. Data were 
reduced with the HEAsoft package v6.14 and the 20130118 
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Figure 5. X-ray light curve in the 2-10 keV energy range during the period 
1999-2012. Square: Chandra', triangle: XMM-Newton', circles: Swift-XKT. 

CALDB-UVOTA release. We extracted the source counts from 
a circle with 5 arcsec radius centred on the source. Background 
counts were derived from an annular region centred on the source 
with 15 and 20 arcsec radii. The observed magnitudes are reported 
in Table|3 During 2008-2011 the difference between the maximum 
and minimum magnitude is 0.7, 0.9, 0.7, 1.3, 1.3 and 1.4 mag from 
the V to the uvw2 band, with the peak of the activity observed on 
2010 June 19. 

As in e.g. iRaiteri et alJ ( 1201 ih . we calculated the effective 
wavelengths, count-to-flux conversion factors (CFa), and amount 
of Galactic extinction in the UVOT bands (Aa) by convolving 
the physical quantities with a power-law fit to the source flux and 
with the UVOT effective areas. In particular, to derive t he Aa we 
adopted the mean extinction law bv ICardelli et al.l d 19891) with Rv 
= AvlE{B — U) = 3.1 (the standard valu e for the diffuse inter stel- 
lar medium), as well as Av = 0.112 from ISchlaflv et all ) l201lh . To 
obtain the UVOT de-reddened fluxes we multiplied the count rates 
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Table 3. Results of SwiftAJVOI observations of PKS 0521—36 in observed magnitude. 


Date 

(MJD) 

Date 

(UT) 

V 

(mag) 

b 

(mag) 

u 

(mag) 

wl 

(mag) 

m2 

(mag) 

w2 

(mag) 

53516 

2005-26-05 

- 

- 

- 

15.94±0.02 

- 

- 

54498 

2008-02-02 

- 

- 

- 

15.11±0.02 

- 

- 

54503 

2008-02-07 

- 

- 

- 

15.67±0.04 

- 

- 

54504 

2008-02-08 

15.16±0.03 

15.89±0.03 

15.41±0.03 

15.58±0.03 

15.71±0.03 

15.82±0.03 

54509 

2008-02-13 

15.15±0.03 

15.93±0.03 

15.42±0.03 

15.59±0.03 

15.69±0.04 

15.83±0.03 

55260 

2010-03-05 

14.94±0.04 

15.68±0.03 

15.13±0.03 

15.34±0.03 

15.44±0.04 

15.59±0.03 

55263 

2010-03-08 

- 

- 

- 

15.18±0.02 

- 

- 

55366 

2010-06-19 

14.58±0.03 

15.12±0.02 

14.49±0.03 

14.54±0.02 

14.63±0.03 

14.54±0.03 

55370 

2010-06-23 

14.98±0.03 

15.67±0.03 

15.18±0.03 

15.36±0.03 

15.51±0.03 

15.63±0.03 

55382 

2010-07-05 

14.83±0.03 

15.48±0.03 

14.93±0.03 

15.11±0.03 

15.17±0.03 

15.27±0.03 

55386 

2010-07-09 

14.93±0.03 

15.61±0.03 

15.08±0.03 

15.28±0.03 

15.33±0.03 

15.46±0.03 

55390 

2010-07-13 

15.09±0.03 

15.87±0.03 

15.33±0.03 

15.44±0.03 

15.67±0.04 

15.74±0.03 

55626 

2011-03-06 

15.27±0.05 

16.05±0.04 

15.58±0.04 

15.71±0.04 

15.90±0.05 

15.97±0.04 

55627 

2011-03-07 

- 

- 

- 

15.74±0.02 

- 

- 

55661 

2011-04-09 

15.14±0.05 

15.95±0.04 

15.50±0.05 

15.59±0.05 

15.76±0.06 

15.83±0.04 

55662 

2011-04-10 

15.10±0.05 

15.90±0.04 

15.39±0.04 

15.56±0.05 

15.69±0.05 

15.83±0.04 

55804 

2011-08-31 

- 

- 

15.63±0.03 

- 

- 

- 

55805 

2011-09-01 

15.23±0.05 

15.99±0.04 

15.53±0.04 

15.81±0.05 

15.70±0.05 

15.85±0.03 

55809 

2011-09-05 

14.94±0.04 

15.57±0.03 

15.13±0.04 

15.39±0.04 

15.58±0.04 

15.78±0.03 


by the CFa and corrected for the corresponding Galactic extinction 
values Aa. 

3.2 XMM-Newton and Chandra observations 

PKS 0521-36 was observed by the Chandra and XMM-Newton 
X-ray satellites on 1999 December 31 (MJD 51543) and 2002 
May 26 (MJD 52420), respectively. The spectral analysis was 
performed using the XSPEC vl2.7 package. 

Chandra ACIS-S spectra and instrument responses were gen¬ 
erated using the Chandra Interactive Analysis of Observations 
(CIAO) v4.3 and the related Chandra calibration database. The 
0.3-7 keV image shows the presence of an unresolved core, a dif¬ 
fuse ha lo and a jet feature coi ncident with the inner radio/optical 
jet (see lBirkinshaw et alj2002ll . 

The nuclear spectrum was extracted from an annular region 
(Pn = 0.3 arcsec, lout = 1.2 arcsec) in order to minimize pile up 
effectfl The background counts were extracted from three source- 
free circular regions of r = 1.2 arcsec each. The 0.5-7 keV spectrum 
is well fitted (x^/dof=187/199) by a power law absorbed by a col¬ 
umn density slightly in excess of the Galactic value, Nh = (4.82 ± 
0.03)xl0^° cm"^ 

The spectrum of the jet was extracted from a circle of radius 
1 arcsec at a distance of 1.2 arcsec from the core. The background 
was chosen from three source-free circular regions. The best fit 
model is a power law with photon index Tx = 2.2 ± 0.3. The 2- 
10 keV flux of the jet is roughly twenty times lower than the core 
flux. Chandra spectral results for the nucleus and the jet are listed 
in Table m 

The XMM-Newton EPIC observation of PKS 0521-36 was 
analyzed using the Science Analysis System (SAS) vl 1.0 software 
and available calibration files. Time intervals affected by high 
background were excluded. After this data cleaning, we obtained a 

® Using the PIMMS software I http://cxc.harvard.edu/toolkit/pimms.jsp I 
the pile up for this source was estimated to be around 4 per cent 


net exposure of 26.6 ks and a count rate of 3.50 ± 0.01 count s“ 
for the pn. The source and background spectra were extracted from 
circular regions of 32 arcsec and 44 arcsec radius, respectively. 
The response matrices were created using the SAS commands 
RMFGEN and ARFGEN. The nuclear data are not piled up. Data 
were grouped into 25 counts per bin in order to apply the 
statistic. A power law, absorbed by Galactic column density, plus 
a thermal (APEC) co mponent is a good para metrization of the data 
(X^/dof=972/925). In iFoschini et al.l (l2006l) . a broken power law 
model absorbed by the Galactic column was reported as the best 
fit to these XMM-Newton data. However the 0.5-10 keV best fit 
model turns out to be a combination of a broken power-law, de¬ 
scribing the nuclear radiation, plus a Raymond-Smith component 
(X^/dof=947/923), which accounts for the thermal emission of the 
diffuse halo around the source, necessarily inside the 32 arcsec 
extraction radius. It is interesting to note that the two spectral 
indices of the broken power-law are similar to those observed for 
the nucleus and the jet by Chandra. This may indicate that we are 
observing both these components in the XMM-Newton spectrum. 
An alternative possibility is that in the XMM-Newton spectrum we 
are observing the tail of the synchrotron emission and the rise of 
the IC component from the jet. There is no significant evidence for 
a Fe Kq emission line. XMM-Newton spectral results are listed in 
Tabled 

4 RADIO DATA 
4.1 VLB A observations 

PKS 0521—36 was observed with the VLBA at 8 and 15 GHz on 
2003 January 7. The central frequencies of the two bands were 
8.421 GHz and 15.356 GHz and the total bandwidth was 32 MHz 
for each band. The on-source time was about 30 and 50 min at 8 
and 15 GHz, respectively, divided into 10 scans distributed over a 
wide range of hour angle, to ensure the best sampling of the (u, v)- 
plane. One station (North Liberty) could not provide data, while 
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Table 4. Spectral results of the nucleus and the jet of PKS 0521—36 for the Chandra observation performed on 1999 December 31. A power-law model was 
used. 


Component 

Parameter 

Value 

Nucleus 

Nh (cm-2) 

(4.82 ± 0.03) xlO^o 


r 

1.7 ±0.1 


norm (ph keV”^ cm“^ s~^) 

(3.1 ±0.4)x10-4 


F2-l0keV (erg cm-2 g-l) 

(1.3 ±0.2)x 10-12 

Jet 

Nh (cm-2) 

3.58x1020 (fixed) 


r 

2.2 ± 0.3 


norm (ph keV”^ cm~^ s~^) 

(2.4tJ'2)xlO-5 


F2-l0keV (ergcm-2 g-l) 

(6.0t?;g)x 10-11 


Table 6. Model-fit parameters. Columns 1 and 2: flux density S of the components; columns 3 and 4: the radial distance d from the core at 8 and 15 GHz, 
respectively; columns 5 and 6: the polar angle 9, measured north to east at 8 GHz and 15 GHz, respectively; columns 7 and 8: the major axis a of the component 
(0 indicates a delta component); column 9: the ratio b/a between the minor and major axis of the component (1 indicates a circular Gaussian); column 10: the 
orientation cf) of the major axis on the plane of the sky; column 11: the spectral index ar. Numbers in italics refer to parameters held fixed in the model fit. 
“These values are left free to vary in the model fit at 8 GHz and fixed at 15 GHz. 


Ss 

(mJy) 

Sl5 

(mJy) 

ds 

(mas) 

^^15 

(mas) 

ds 

(°) 

016 

(°) 

as 

(mas) 

ai5 

(mas) 

hfa 

4> 

(°) 

Qr 

545 

815 

0.00 

0.00 

0.0 

0.0 

0.00 

0.00 

0 

0 

- 0.7 

414 

584 

0.28 

0.27 

-54.9 

-33.6 

0.37 

0.27 

1 

0 

- 0.6 

233 

195 

0.96 

0.96 

-42.9 

-43.0 

0.44 

0.25 

1 

0 

0.3 

100 

72 

2.06 

2.28 

-34.1 

-42.8 

0.87 

1.11 

1 

0 

0.5 

33 

12 

4.30 

4.35 

-40.8 

-39.9 

1.18 

1.18 

1 

0 

1.7 

38 

25 

10.8 

10.8 

-55.5 

-55.5 

4.58 

4.58 

1 

0 

0.7 

298 

220 

27.9 

27.9 

-45.9 

-45.9 

7.95 

7.95 

0.58“ 

-19.7“ 

0.5 

94 

59 

35.3 

35.3 

-47.3 

-47.3 

8.53 

8.53 

0.47“ 

-6.9“ 

0.8 


Table 5. Spectral results for the XMM—Newton observation of PKS 
0521—36 performed on 2002 May 26. A broken power-law model with 
A^h fixed to the Galactic absorption plus a Raymond-Smith component was 
used. 


Parameter 

Value 

Ti 

2.4 ± 0.2 

Eb (keV) 

G. /o_o.04 

r2 

1.7 ± 0.2 

norm (phkeV“^ cm~^ s~^) 

(1.7 ± 0.1)xlO-o 

kT (eV) 

U.o3_0 07 

norm (phkeV^ cm~^ s“^) 

( 8.2 ±2.7)x 10-5 

F 2-10 keV (erg cm-2 s-l) 

(8.0 ±0.7)x 10-12 


two performed poorly because of snow (Hancock) and hardware 
problems (Brewster). 

We downloaded the data from the National Radio Astronomy 
Observatory (NRAO) archive and performed the typical steps of 
calibration in the astronomical image processing system (AlPS). 
We used system temperature and gain curve tables for the visibility 
amplitude calibration and FRING for the delay, rate, and visibility 
phase. Uncertainties on the amplitude calibration are about 10%. 
Several cycles of self-calibration were then performed in DIFMAP 
to improve the initial calibration, both in phase and amplitude, and 
produce final images. 

In Figs. [6] and [7] we show total intensity images of 
PKS 0521—36 at 8 and 15 GHz, respectively. The source is well 


detected at both frequencies, with a bright compact core and a one¬ 
sided jet oriented in PA —45° on the plane of the sky. The bright¬ 

ness profile along the jet axis decreases rapidly with increasing dis¬ 
tance from the core, but it then suddenly rises again at ~ 30 mas. 
The total flux density is 1.78 and 1.93 Jy at 8 and 15 GHz, respec¬ 
tively, with a peak brightness of 0.94 and 1.21 Jy beam~ ^. The core 
has an inverted spectral index of = -0.4±0.2 (S(v) oc 
The bright jet emission region at ~ 30 mas has a flux density 
Ss ~ 405 mjy and S 15 ~ 265 mjy, corresponding to a spectral 
index of a = 0.7 ± 0.2. 


In Table|^ we report the parameters of a model fit with Gaus¬ 
sian components to the visibility plane. We have first model-fitted 
the 8 GHz visibilities with delta (for the core), circular (for the inner 
jet), and elliptical (for the outer blob) Gaussian components, until 
we reached convergence. Then, we used this model as a starting 
condition for the model fit to the 15 GHz visibilities. In this lat¬ 
ter fit, we kept the position and size of the outer three components 
fixed, while all the other parameters (including the flux density of 
every component) were left free to vary. 

In order to maximize sensitivity to the diffuse jet emission, 
we also produced an image applying a Gaussian taper of 0.3 at a 
radius of 30 MA to the 8 GHz dataset. This yields a resolution of 
4.9 mas x 14.2 mas in PA 2?3, a peak brightness in the jet region of 
about 185 mJy beam~^, and an off source peak of 1.4 mJy beam“^ 
on the counter-jet side. 
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Figure 6. VLBA image of PKS 0521—36 at 8 GHz. On the image we pro¬ 
vide the restoring beam, plotted in the bottom-left corner, the peak flux 
density in Jy/beam, and the first contour (f.c.) intensity in mJy/beam, which 
is three times the off-source noise level. Contour levels increase by a factor 
of 2. 

VLBA 15 GHz 

peak =1.21 Jy/beam; f.c. = 2.4 mJy/beam 



I 

20 10 0 -10 -20 
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Figure 7. VLBA image of PKS 0521—36 at 15 GHz. On the image we 
provide the restoring beam, plotted in the bottom-left comer, the peak flux 
density in Jy/beam, and the first contour (f.c.) intensity in mJy/beam, which 
is three times the off-source noise level. Contour levels increase by a factor 
of 2. 


4.2 VLA data 

To obtain information on the radio flux density on kiloparsec-scale, 
we retrieved archival NRAO VLA data at 4.9 and 8.4 GHz. Ob¬ 
servations of PKS 0521—36 were performed at 4.9 GHz on 2009 
August 6, when the array was in C configuration, while 8.4-GHz 
observations were carried out on 2009 October 3 when the array 
was in the hybrid C and D configuration. During both runs the 
target was observed for about 40 s. A few antennas did not par¬ 


Table 7. VLA flux density and spectral index of PKS 0521—36. 


Comp 

54.9 

S8.4 

^8.4 

^4.9 


Jy 

Jy 


c 

3.72±0.18 

3.94±0.20 

-0.1 ±0.1 

s 

2.18±0.11 

1.35±0.07 

1.0±0.1 

Ext 

1.12±0.06 

0.35±0.03 

2.1±0.2 

Tot 

7.02±0.35 

5.64±0.28 

0.4±0.1 


ticipate in the observations due to a system upgrade. The primary 
flux density calibrator was 3C48. Uncertainties on the amplitude 
calibration (la) are about 5 per cent at both frequencies. The data 
reduction was carried out following the standard procedures for the 
VLA implemented in the NRAO AlPS package. The flux density 
at each frequency was measured on the final image produced after 
a few phase-only self-calibration iterations using TVSTAT, which 
performs an aperture integration over a selected region on the image 
plane. Flux density of the compact components was measured by 
the AlPS task JMFIT, which performs a Gaussian fit of the source 
component on the image plane. Flux densities are reported in Ta- 
ble|3 The resolution achieved is about 17 arcsec x 4 arcsec and 11 
arcsec x6 arcsec at 4.9 and 8.4 GHz, respectively. 

In Figs. and we show the total intensity images at 4.9 
and 8.4 GHz, respectively. The radio emission is dominated by 
two compact components, labelled C and S, whose flux density is 
about 80 per cent and 73 per cent of the source total flux density 
at 4.9 and 8.4 GHz, respectively. The central component, coinci¬ 
dent with the optical nucleus, has an inverted spectrum (alig ~ 
-O.liO.l) and hosts the source core. This value is larger than that 
derived from VLBA data and this is likely due to the contribution 
of the jet that cannot be separated from the core component due to 
the lower resolution of the VLA observations. A jet-like structure, 
labelled J, emerges from the core with a position angle of about 
—45°, in agreement with what is found in high-resolution VLBA 
images (Fig. IT}. At about 5 arcsec (5 kpc) from the core, the jet 
bends toward the west, producing an extended structure of about 
30 arcsec (~ 30 k pc) in size, in agreement with what i s found at 
lower frequencies dFalomo et al.ll200^ : lEkers et al.lll989h . Compo¬ 
nent S is located about 7 arcsec (~ 7 kpc) from the central compo¬ 
nent with a position angle of 120° and it marks the location of the 
hotspot. Its spectral index is rather steep a® I = l.OiO.l and it is 
likely due to the low resolution of these observations, which does 
not allow us to disentangle the contribution of the hotspot from that 
of the lobe. It is worth noting that the spectral indices have been 
computed assuming the total flux density measured on images with 
different resolution, which may produce an artificial steepening of 
the spectrum. The errors on the spectral index have been computed 
by means of the error propagation theory. 

4.3 UMIfAO data 

UMRAO centimeter band total flux density observations were 
obtained with the University of Michigan 26 m paraboloid lo¬ 
cated in Dexter, Michigan, USA. The instrument is equipped with 
transistor-based radiometers operating at frequencies centered at 
4.8, 8.0, and 14.5 GHz with bandwidths of 0.68, 0.79, and 1.68 
GHz, respectively. Dual horn feed systems are used at 8 and 14.5 
GHz. Each observation consisted of a series of 8 to 16 individual 
measurements over approximately a 25 to 45 min time period, uti- 
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VLA 4.9 GHz 



Figure 8 . VLA image of PKS 0521—36 at 4.9 GHz. On the image we pro¬ 
vide the restoring beam, plotted in the bottom-left comer, the peak flux 
density in Jy/beam, and the first contour (f.c.) intensity in mJy/beam, which 
is three times the off-source noise level. Contour levels increase by a factor 
of 2. 

lizing an ON-ON technique (switching the target source between 
the two feed horns which are closely spaced on the sky) at 8.0 
and 14.5 GHz. As part of the observing procedure, drift scans were 
made across strong sources to verify the telescope pointing correc¬ 
tion curves, and observations of nearby calibrators were obtained 
every 1 to 2 hours to correct for temporal changes in the antenna 
aperture efficienc y. Further details about UMRAO are reported in 
lAller et alj j2014|) . UMRAO data at 8 GHz and 14.5 GHz are rep¬ 
resented in Fig.llOl 


4.4 SMA data 

The 230 GHz (1.3 mm) light curve was obtained at the SMA on 
Mauna Kea (Hawaii). PKS 0521—36 is abright AGN included in an 
ongoing monitoring program at the SMA to determine the fluxes of 
compact extragalactic radio sources that can be used as calibrators 
at mm waveleng ths. Details of t he ob servations and data reduction 
can be found in iGurwell et alj ( l2007h . Data from this programme 
are updated regularly and are available at the SMA websitqj SMA 
data at 230 GHz are shown in Fig. 


5 DISCUSSION 

5.1 Identification of the 7 -ray source 

The association of a 7 -ray source with a low-energy counterpart 
is fundamental for understanding its physical properties. The ac¬ 
curacy with which most of the 7 -ray sources were localized with 

^ http://smal.sma.hawaii.edu/callist/callist.html. Use in publication re¬ 
quires obtaining permission in advance. 



Figure 9. VLA image of PKS 0521—36 at 8.4 GHz. On the image we pro¬ 
vide the restoring beam, plotted in the bottom-left comer, the peak flux 
density in Jy/beam, and the first contour (f.c.) intensity in mJy/beam, which 
is three times the off-source noise level. Contour levels increase by a factor 
of 2. 


EGRET was not enough to associate them with a known counter¬ 
part, leaving the legacy of a large fraction of unidentified sources in 
7 -rays. Thanks to larger photon statistics and unprecedented angu¬ 
lar resolution at high energies Fermi-LAT is able to investigate the 
natures and counterparts of these unidentified sources, although the 
7 -ray activity state as well as correlated flux variability at different 
frequencies play an important role for the ide ntification process of 
a 7 -ray source (e.g.. lD’Ammando et al.ll2oT2h . 

A 7 -ray source at a position compatible with PKS 0521—36 
was detected with a 4.3a significance by EGRET during 1992 May 
14—June 4, with a flux of (1.8 ± 0.5)xl0~^ p h cm“^ s“^ and a 
photon index F-, = 2.16 ± 0.36 jLin et al.ll99A . As a result of that 
activity the source was included in the Seco nd EGRET catalogue as 
2EG J0524-3630 dThompson et al.lll995h . The 7 -ray source 3EG 
J0530—3626 was included in the Third EGRET catalogue, but PKS 
0521—36 was outside the 99 per cent confidence contour in that 
catalogue, and the associ ation between the 7 -rav source and the 
AGN was not confirmed jHartman et al.|[l999l) . The same conclu¬ 
sion was drawn for the 7 -ray source EGR J0 529—3608 reported in 
the rev ised EGRET catalogue proposed by ICasandiian & GrenieJ 

(l2008h . 

No 7 -ray source associated with PKS 0521—36 was reported 
in the LAT Bright S ource List obtained after the first three months 
of Fermi operation jAbdo et^l2009ll . but the AGN was associ¬ 
ated with lEGL J0522.8-3632, 2PGL J0523.0-3628 and 3PGL 
J0522.9—3628 in the following Fermi c atalogues jAbdo et^ 
l2010al:lNolan et al.ll2012l:lAcero et al.ll2015h . 

Most of the associations of 7 -ray sources with AGN in the 
Fermi catalogues have high probability of being true, but a firm 
identification needs to be confirmed by means of correlated vari¬ 
ability in different energy bands. The contemporaneous increase of 
activity observed in the optical, UV, X-ray, and 7 -ray bands in 2010 
June (Pig. [M is a distinctive signature for the identification of the 
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Figure 10. Multifrequency light curve for PKS 0521—36. The period covered is 2008 August 4-2012 August 4 (MJD 54682-56143). The data were collected 
(from top to bottom) by Fermr-LAT ( 7 -rays; in units of I0“® ph cm~^ s“^), Swift-XKT (X-rays; in units of 10“^^ erg cm~^ s~^), Swift-WOT (tul, m2, 
and w2 filters; in mJy), (v, b, and u filters; in mJy), SMA (230 GHz; in Jy), UMRAO (8 and 15 GHz, triangles and squares, respectively; in Jy). The LAT data 
were obtained with the photon index left free to vary (see Fig.l^. 


7 -ray source with PKS 0521—36. In the following, we investigate 
the properties of PKS 0521—36 from radio to 7 -ray energy bands 
with the goal of unveiling the nature of this 7 -ray emitting AGN. 


5.2 Gamma-ray properties 

By comparing the 7 -ray properties of PKS 0521—36 with those 
of the different types of AGN detected by LAT, we noted that the 


average photon index of this source (T^ = 2.45 ± 0.02) is simi- 
lar to the average valu e observed for FSRQ (T-, = 2.44 ± 0.20; 
Ackermaim et^ r20 1 5h and misaligned AGN (T-, = 2.42 ± 0.28; 


Abdo et all[ 2 ^ . As shown n Fig. there is no general trend 


between 7 -ray flux and photon index for PKS 0521—36, but a 
significant hardening of the spectrum was observed during the 7 - 
ray flaring activity in 2010 June. This behaviour is in common 
with bright FSRQ and low-synchrotron-peaked BL Lacs detected 
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by LAT dAbdo et alj 12010 j) . The hardening of the spectrum to¬ 
gether with the low redshift make this source a promising target 
for the current (i.e., MAGIC, VERITAS, and H.E.S.S.) and next 
generation of imaging atmospheric Cherenkov telescopes (e.g., the 
Cherenkov Telescope Array). It is interesting to note that one of the 
few radio galaxies detected at Very Hi gh Energy (VHE; E > 100 
GeV), NGC 1275 dAleksic et al.ll201^ . showed flux and spectral 
index similar to those of PKS 0521—36 during the 2010 June flare. 

Contrary to the photon index, the average apparent isotropic 
7 -ray luminosity of PKS 0521—36 (L-^ = 5.2x 10^^ erg s“^) seems 
to lie in the region usual ly occupied by BE Lacs and radio galaxies 
dAckermann et al.ll2015l) . However, the lower apparent luminosity 
of PKS 0521 —36 with respect to both FSRQ and the steep spectrum 
radio quasars (SSRQ) 3C 380 and 3C 207 may be due to the close 
proximity of the source. In fact, a 7 -ray flaring activity with peak 
flux of 1-2x10“® ph cm“^ s“^, as observed for PKS 0521—36, 
is quite uncommon in BL Lacs (e.g.. ICannon et ai] l 2010 h and in 
misaligne d AGN, the on ly exception being the radio galaxy NGC 
1275 (e.g.. lCiprinill2013h . 

It is not clear whether the most powerful jets of the two parent 
populations (FSRQ and FRII radio galaxies) have analogous struc¬ 
tures. The difference in the Jet structure may be related to the differ¬ 
ent environment or to jet properties. We noted that the FR II radio 
galaxy Cygnus A, one of the most powerful radio sources with a 
redshift comparable to that of PKS 0521—36, has not been detected 
in 7 -rays so far. On the other hand, Pictor A, an FR II radio galaxy 
at red shift z = 0.035, was detected by Fermi-hAI jAcero et al.l 
l2015t) . but with an apparent 7 -ray luminosity that is two order of 
magnitudes lower than that observed for PKS 0521—36. This may 
support the scenario of a complex jet structure in radio galaxies 
(e.g., velocity gradients, spine-layer) with possible differences in 
the emission mechanisms at higher energies also among objects of 
the same class. Alternatively, the different viewing angle with re¬ 
spect to the observer, and thus the different Doppler boosting of the 
jet emission, may lead to the detection or not of radio galaxies in 
7 -rays. 


No significant difference was observed between the average 
7 -ray flux detected by EGRET in 1992 and LAT during 2008- 
2012. This is different from the long-term variability observed 
for the radio-galaxy NGC 1275 between the EGRET and Fermi 
era. In that case, the early LAT observations revealed a flux at 
least seven times hig her than the upper limit obtained by EGRET 
dKataoka et alKOlOl) . Of particular interest is the short variability 
observed during flaring periods with a doubling time scale of the 
order of a few hours, co mpatible with the variability observed in 
some brigh t FSRQ (e.g.. ISaito et al.|[2ol3l : iTavecchio et al ]H; 
lAbdo et ^l2010cl) . If the emitting region fills the entire cross sec¬ 
tion of the jet, this rapid variability suggests a compact 7 -ray emit¬ 
ting region and the location of the 7 -ray emission may be within 
the broad-line region (BLR). Following the causality argument and 
considering the 12 -hr variability observed during the 7 -ray flares 
the size of the emitting region should be i? < cf„ar5/(l + z) 
= 6.3 xlO^® cm (assuming <5 = 5; see Section [531 i. This size is 
smaller than the values usually inferred from the SED modeling 
of blazars w ithin the framework o f the standard one-zone leptonic 
model (e.g.. iGhisellini et'H]|201 ih . Alternatively, episodes of fast 
varia bility can be produced in a jets-within-iet/m ini-jets scenario 
(e.g. IGhisellini et S]|2008l : iGiannios et alj|2009l) and in this case 
the location may b e further out, as proposed e.g. for M87 (e.g. 
iGiroletti et al.l20 1^. or pro duced by turbulent cells in the relativis¬ 
tic plasma ( Marscheill2014h . 


5.3 Radio properties 


On a kpc scale, the radio source shows a two-sided structure where 
jets, lobes, and a hot spot are clearly visible in addition to the 
core. The structure revealed by the V LBA images agrees with 
those reported in p revious works (e.g. iTingav & Edwards! 12002 : 
IGiroletti et al.ll20()4) . However, since the data were taken at distant 
epochs and at different frequencies, it is difficult to make a reli¬ 
able identification of the individual components. The best possible 
match is probably on the bright diffuse region at ~ 30 mas, whose 
flux density and position are consistent in the various datasets, 
and essentially stationary within the uncertainty related to its large 
size (~ 8 mas). This bright resolved structure is reminiscent of 
two other 7 -ray emitti ng AGN somewh at different fro m the typi- 
cal blazars, i.e. M87 1 Abdo et^l2009h and 3C 120 jAbdo et^ 
l201(llKataoka et alj|201 ll) . On parsec scales, both sources show a 
one-sided jet whose brightness profile along the axis initially de¬ 
creases down to the noise level and then steeply rises in an ex¬ 
tended, stationary feat ure. The projected distan ce of this feature 
is ~ 50 pc in 3C 12 0 1 Roca-So gorb^t^ 2010l) and ~ 70 pc in 
M87 dGiroletti et alj|2012l:lcheung et alj|200’^ ln PKS 0521-36, 
a bright, extended jet feature is found at nearly 30 pc (projected). 
These features may be interpreted as recollimation shocks, and at 
least in the case of M87 it has been proposed that high-energy 
emission up to the TeV regi me might originate at this location 
dBromberg & LevinsonI l2009l) . Recently, a similar scenario was 
sugges ted also for the rapid V HE flaring activity of the FSRQ 4C 
21.35 Iayecchio_gt_alJ| 201 lJ), although different models are pro¬ 


posed I Ackermann et alj2014l : lDermer et al.l2012l : lTavecchio et al.l 

|2012h . 

We estimated the ranges of viewing angles 9 and jet velocity 
P (= v/c) from the jet/counter-jet brightness ratio from the radio 
images presented in Sec. l4.ll Assuming the off-nuclear peak on the 
counter-jet side as an upper limit to any possible counter-jet bright¬ 
ness, we calculate a jet-counterjet brightness ratio R ~ 130. In the 
standard hypothesis that jet and counter-jet are intrinsically similar 
and that any asymmetry is entirely due to relativistic beaming, this 
places limits on the jet velocity and orientation by 


^ . 1 + Pcos9 .2+cc 
Ecj t — P cos6 


( 1 ) 


and we obtain P > 0.75, <5 > 1.5, 9 < 41°, assuming a = 0.7. 

The radio spectrum below 1 GHz is dominated by the flux den¬ 
sity arising from the exte nded struct ures. The spectral index com¬ 
puted between 160 MHz 1 Sled 197% a nd 1.4 GHz from the NRAO 
VLA Sky Survey jCondon et al.lll99d) results to be 0.7. This radio 
spectral shape and the optical emission suggest that PKS 0521—36 
may be an intermediate object between a broad-line radio galaxy 
(BLRG) and an SSRQ. 

We computed the core dominance (CD) of PKS 0521—36 follow¬ 
ing the definition in lOrr & Brownel ( ll982l) : lAbdo et al.l ( l2010bl) : 


CD 


5c, 


5tot 5c, 


( 2 ) 


where the core flux density Score and the total flux density 5tot 
refer to the source rest frame (5 = 5obs x For the core 

flux density we consider the value measured on the VLA image at 
4.9 GHz and reported in Table 6 , and we assume a flat spectrum 
a = 0. In order to prevent missing flux density associated with 
extended emission, we consider the total flux density measured by 
the Parkes telescope at 5 GHz, which is 5tot = 8.1 Jy ll Wright et al.l 
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Il996h . and we assume a spectral index a = 0.7. We obtain log(CD) 
= —0.10. This value is similar to those shown by the 7 -rav emittin g 
SSRQ 3C 380 and 3C 207 and BLRG 3C 111 jAbdo et al.ll2010bh . 

In order to estimate the radio power and to make a comparison 
with the 7 -ray emitting misaligned AGN, we computed the radio 
luminosity at 178 MHz by using: 

L = 4TvdlS{l + (3) 

where is the luminosity distance, and S is the flux density at 
178 MHz, which has been extrapolated from the values reported in 
the literature and corresponds to 67.5 Jy (source rest frame). The 
corresponding luminosity is L = 4.8 x 10^® W/Hz, which is a few 
orders of magnitude lower than the luminosity of the SSRQ 3C 380 
and 3C 207, but is comparable with the luminosity of the BLRG 3C 
111 . 

5.4 X-ray to mm behaviour 

During 2005-2011 the source was monitored by the Swift satellite 
only sporadically, but a quite prominent flare was observed from 
optical to X-rays in 2010 June, contemporaneous with the peak of 
7 -ray flaring activity (see Section ED. The small variability am¬ 
plitude (calculated as the ratio of maximum to minimum flux) ob¬ 
served in X-rays (-^ 3) with respect to that observed in 7 -rays (~ 13 
on a monthly time-scale up to a factor of ~ 50 on a daily time-scale) 
together with the lack of hardening of the spectrum may be an In¬ 
dication that the X-ray emission is produced by the low-energy tail 
of the same electron distribution, as o bserved in FSRQ (e.g., PKS 
1510-089: [P’Ammando etaD|201lh . On a time-scale of years a 
change by a factor of ~10 in flux was detected between the Chan¬ 
dra and Swift observations. A second less prominent increase of the 
X-ray activity was observed on 2011 September 5, with no evident 
counterpart in 7 -rays. This may be due to different mechanisms or 
emission sites that give significant contribution in X-rays and 7 - 
rays in that period. The photon indices observed by Swift-BAT and 
Swi/r-XRT in hard X-rays (T x ~1.8) and soft X-rays (T x -^1.7) 
indicates that the peak of the IC emission should be at MeV ener¬ 
gies, according to the spectra observed by Termi-LAT (T-, > 2 ). 

The variability amplitude in the optical-UV bands covered by 
Swift-VWOT is V ~1.9, b -2.4, u -2.9, wl -3.2, m2 -3.2, and 
w2 —3.7. We observed a larger variability amplitude at higher fre¬ 
quencies, contrary to what is observed in quasars, where the contri¬ 
bution from the thermal accretion disc significantly in creases in the 
UV p art of the spectrum, diluting the jet emission (e.g. lRaiteri et al.l 
l2012h . This suggests the lack of a prominent accretion disc in the 
broad-band spectrum of the source. We noted that during the flaring 
activity in 2010 June the optical and UV emission seems to peak a 
few days before the X-ray peak. This may be due the coverage pro¬ 
vided by the Swift observations, and thus the lack of an observation 
at the time of the daily 7 -ray peak, or to the complex jet structure 
with different parts of the jet responsible for the emission at differ¬ 
ent frequencies. 

The SMA data at 230 GHz show an amplitude variability 
of —4, with a peak flux density of 8 Jy on 2010 December 30 
(MID 55560). This event at millimeter wavelengths may be the de¬ 
layed flare with respect to the 7 -ray flare observed in 2010 June. 
The increase of flux density was observed also by UMRAO at 15 
and 8 GHz, with the difference in time of about 1 month and 2 
months, respectively, with respect to the 230 GHz emission likely 
related to opacity effects due to synchrotron self-absorption. Re¬ 
cently, [Fuhmiann_et_^ ( 120 14h investigated the correlation between 


radio and 7 -ray emission for a sample of bright Fermi blazars. Per¬ 
forming a stacking analysis of the objects in the sample, they found 
at 142 GHz a time-lag of (7±9) days between the radio and 7 -ray 
emission, with the 7 -ray leading the radio emission. In this context 
a delay of six months between the peak at 230 GHz and the 7 -ray 
one for PKS 0521—36 is extremely large and thus unlikely, open¬ 
ing alternative scenarios. The emitting region responsible for the 
optical-to- 7 -ray emission may be different from the region produc¬ 
ing the mm and radio emission. However, we note that SMA data 
are not available in 2010 May-June, and therefore we cannot rule 
out a mm flare in the same period of the 7 -ray flare. 

5.5 SED modeling 

We build two SED for PKS 0521—36 in two different activity 
states: the high state in 2010 June and a low state in 2010 February- 
March. The flaring SED includes the Fermi-LAT spectrum built 
with data centred on 2010 June 22-July 5 (MJD 55369-55382), 
the optical, UV, and X-ray data collected by Swift on 2010 June 23 
(MJD 55370), and 230 GHz data collected by SMA on 2010 July 1 
(MJD 55378). The low SED includes the EAT spectrum built with 
data centred on 2010 Eebruary 1-March 31 (MJD 55228-55286), 
the optical, UV, and X-ray data collected by Swift on 2010 March 5 
(MJD 55260), and 230 GHz data collected by SMA on 2010 Eebru¬ 
ary 24 (MJD 55251). 

The modeling of the SED of PKS 0521—36 within the frame¬ 
work of the standard one-zone leptonic model provides results 
at odds with the bulk of the blazar population. In particular, a 
quite low Doppler fac tor is required (5 ~ 3; e.g., IPian et aOl 19961 : 
iGhisellini et al.ll201lh . implying a relatively large viewing angle 

^ 15° and a low bulk Lorentz factor. This result can be 
traced back to the relatively small separation between the two 
bumps in the SED. This conclusion is corroborated by several hints, 
such as the existence of a large-scale optical jet and the moder¬ 
ate core-to-extended radio flux ratio. In the case of a relatively 
large viewing angle and thus small beaming amplification, it is 
not excluded that also slow regions of the jets contribute to the 
observed emission. This would be the case in the framework of 
the structured jet scenario, in which the flow is supposed to be 
composed of a fast inner core - or spine - surrounded by slower 
plasm a - the layer dGhisellini et al.lE005l : iTavecchio & Ghisellinil 
l2008h . Indeed such a scheme has been applied to interpret the emis¬ 
sion of 7 -ray emitting radio galaxies, for which the emission can 
be interpreted as a mix of sp ine and layer components (see also 
ITavecchio & Ghisellinill20i4) . With this motivation, in the follow¬ 
ing we focus our attention on the application of the structured jet 
model to the emission of PKS 0521—36. 

We briefly recall the main features and the parameters of the 
model (refer to Ghiselllni et al. 2005, Tavecchio & Ghisellini 2008, 
2014 for a more complete description). The spine is assumed to 
be a cylinder with height Ha and radius Ra (the subscripts ‘s’ and 
‘F stand for spine and layer, respectively). The layer is modeled as 
a hollow cylinder with height Hi, internal radius Ra and external 
radius Ri = 1.2x J?s. We assume Ha <C Hi, ideally corresponding 
to the case of a perturbation traveling down the spine, surrounded 
by a relatively long and stationary layer of slow plasma. 

In both zones we assume a phenomenological and stationary 
electron energy distribution described as a broken power law with 
the following parameters: minimum, maximum and break Lorentz 
factors 7 min, 7 max and 7 b and indices ni and n 2 . The normaliza¬ 
tion of the electron distribution is parametrized by the emitted syn¬ 
chrotron luminosity, Lsyn- The emitting regions contain a tangled 


























Unveiling the nature of the 'y-ray emitting AGN PKS 0521—36 13 



R 

H 

Tgyn 

B 

7min 

7b 

'7max 

ni 

712 

r 



cm 

cm 

erg s”^ 
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deg 

Layer low 

3.5 X IQi® 

2 X IQi® 

6 X 10^9 

0.5 

100 

3 X 10° 

7 X lO* 

1.5 

3 

3 

6 

Spine low 

3.5 X IQi® 


1.9 X lO"*! 

0.18 

60 

1.2 X 10"* 

2.55 X 10"* 

2 

3.7 

12 

6 

Laver high 

3.5 X IQi® 

2 X IQl® 

7 X 10"*° 

0.5 

100 

3 X 10° 

7 X 10"* 

1.5 

3 

3 

6 

Spine high 

3.5 X IQi® 


2.2 X 10"*! 

0.17 

60 

6 X 10° 

3 X 10"* 

2.1 

3 

12 

6 

Layer low 

5 X IQi® 

2 X IQl® 

1.3 X lO"*! 

1 

80 

3 X 10"* 

10® 

2.75 

4 

3 

15 

Layer high 

5 X IQi® 

2 X IQl® 

5.0 X lO"*! 

1 

100 

10^ 

5 X 10® 

2.7 

4 

3 

15 

Spine* 

5 X IQi® 


4 X 10"*° 

15 

10 

2 X 10° 

6 X 10° 

2.1 

3 

10 

15 

Layer high 

5 X 10^® 

2 X 10^® 

2 X 10"*° 

1 

80 

6 X lO* 

8 X 10® 

2.7 

4 

3 

20 

Spine high 

5 X IQi® 


4 X 10"*"* 

15 

10 

300 

8 X 10° 

2.1 

3 

10 

20 


Table 8. Input parameters of the models for the layer and the spine shown in Fig. \m All quantities (except the bulk Lorentz factors F and the viewing angle 
6v) are measured in the rest frame of the emitting plasma. The external radius of the layer is fixed to the value R 2 = 1.2 x R. *For the case 6v = 15° the 
same spine emission has been used for both states. 


magnetic field Bs, Bi. The relativistic beaming is specified by the 
two Lorentz factors Fs^Ti and by the viewing angle 0^. 

E lectrons emit through s ynchrotron and IC mechanisms. Fol¬ 
lowing [chiseihnret^ | 200 ^, RbUR = 10 ^^ X-^^diL ,45 for 

-L disc — 3x10 ^rg s we obtain -Rblr, — 2x10^® cm. There¬ 
fore the emitting region should be beyond the BLR and for this 
reason the external Compton of the BLR seed photons should be 
negligible. For the IC besides the local synchrotron radiation field 
(synchrotron self-Compton, SSC) we consider the beamed radia¬ 
tion field of the other component. Indeed, because of the relative 
motion, the energy density of radiation produced in one component 
is boosted by the square of the relative Lorentz factor Fret in the rest 
frame of the other, where Frei = rsri(l — /3s/3i). As we will see, 
absorption of 7 -rays through the interaction with the soft radiation 
field, 77 —>■ e''"e“ can be important. In the case under study, radia¬ 
tion is produced and absorbed within the same region and thus the 
“suppression factor” is /(iz)/7o(t^) = (1 —exp[— r-y.y(iz)])/r-,.y(t/) 
which, for large optical depths, is well approximated by 

We note that the number of the parameters is large (almost 
twice that of the simple SSC model). However, the model has to 
satisfy constraints that ca n be used as guidelines in sele cting the 
suitable setup. We refer to lTavecchio & Ghisellinil (l2014h for a de¬ 
tailed discussion. 

Following the indications recalled in the previous subsections 
we consider angles larger than those usually associated with blazar 
jets, < 5°. In Fig. m we show the results of the modeling of 
the SED with three different angles, 9^ = 6 °, 15° and 20°. The 
corresponding parameters are reported in Tabled 

We expect that, for relatively small angles, the spine radiation 
is much more boosted than that of the layer, implying that the spine 
emission component largely dominates the observed SED. This so¬ 
lution is that assumed for the case 6^ = 6 ° (upper panel). The 
relative stability of the low-energy component accompanied by the 
variation which occurred in the IC peak is obtained by changing the 
layer luminosity (and maintaining almost the same parameters for 
the spine). In fact, the spine IC bump is dominated in the high state 
by the component resulting from the IC scattering of the electrons 
in the spine with the photons coming from the layer. Variations of 
the layer luminosity are thus accompanied by substantial variations 
of the IC luminosity of the spine. 

For larger angles the only suitable solution is that the two 
bumps are dominated by the two jet components, i.e. the low en¬ 
ergy bump by the synchrotron component of the spine and the high 


e nergy bump by the IC emissio n of the layer (see the discussion 
in iTavecchio & Ghisellinill2014h . For 9^ = 15° (middle panel) the 
model still reproduces quite well the data. In this case we also as¬ 
sume a stationary spine and a variable layer component. A further 
increase of t he angle (lower panel) is limit ed by a problem already 
discussed in iTavecchio & Ghisellinil l l2014ll . Indeed, to compensate 
the reduced relativistic amplification of the spine emission one has 
to increase its intrinsic luminosity. This radiation field - concen¬ 
trated in the optical-IR band - is the ideal target for the pair pro¬ 
ducing reactions 77 e'^+e~ with 7 -ray photons produced in the 

layer. In Fig.[T2 we show the optical depth and the correspond¬ 
ing “suppression factor” [1 — exp(—r-y.y)]/r-y.y as a function of the 
frequency for the case 9^ = 15° and 9^ = 20°. The shapes of 
the curves, related to the spine target photon spectrum, are similar. 
However, due to the larger intrinsic luminosity, for 9^ = 20° the 
optical depth is larger by a factor of ~ 10. As a consequence, while 
= 15° the source is transparent (r.y-y < 1 up to ~ 20 GeV, cor¬ 
responding to the highest-energy LAT bin in the high state, energy 
marked by the vertical yellow line), for 9^ = 20 °, the optical depth 
reaches unity already at « 1 GeV and rapidly increases, determin¬ 
ing the abrupt cut-off visible in Fig.[TT](lower panel). We therefore 
conclude that, in the framework of the structured jet model, the an¬ 
gle cannot be larger than about 15°. 

In the unification scheme for radio galaxies and blazars we ex¬ 
pect that, once observed at small angles, the SED of a radio galaxy 
resembles that of a typical blazar. Therefore for the remaining two 
cases, = 6 ° and 15° one can further ask what kind of SED 
would the source present if observed at angles more typical for 
blazars. This is shown in Fig. M in which we report the SED as 
recorded by an observer at 9^ — 4°, not changing the other pa¬ 
rameters. For comparison we also display the historical SED of BL 
Lacertae. The case 9^ = 15° (green) results in an SED strongly 
unbalanced (by a factor « 100 ) toward the synchrotron luminosity, 
at odds with the typical blazar SED. This effect is clearly related to 
the fact that the high-energy bump in the case 9^ = 15° is associ¬ 
ated with the layer. At small angles, the layer radiation is much less 
beamed than that of the spine and this determines the substantial 
prominence of the low-energy peak. On the contrary, the SED for 
case 9v = 6 ° (red) resembles that of the prototypical BL Lac. This 
suggests that a viewing angle between 6 ° and 15° is favored. 

In addition in the spine-layer scenario proposed here, the fast vari¬ 
ability observed in 7 -rays is difficult to reconcile with the cases that 
have = 15°, in which the region responsible for the 7 -ray emis- 
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Figure 11. Spectral energy distribution of PKS 0521—36. Red sym¬ 
bols show the high state, green symbols the low state. An'ows repre¬ 
sent 2cr upper limits. Historical data (gray) are from the ASDC archive 
(http://tools.asdc.asi.iti. The three panels report the results 
obtained with the stmctured jet emission model for different values of the 
viewing angle, from top to bottom: 0^ = 6 °, 6v = 15° and 6^ = 20°. The 
dashed and the dotted lines show the emission from the layer and the spine, 
respectively. The thick solid line shows the total. For the case 0v = 15° the 
same spine emission (shown in black) is used for both states. See text for 
details. 

sion is relatively large (5x 10^® cm) and the Doppler factor small 
(d = 4). Reducing the size of the emitting region in these two cases 
should lead to an increase of the compactness of the region and 
thus of the corresponding 7 -ray opacity. This is another indication 
in favour of a relatively small viewing angle. 

However, we cannot rule out the possibility that PKS 
0521—36 is an AGN seen at '^15° and its ‘aligned’ counter¬ 
part is a blazar belonging to a population of synchrotron-dominated 
objects. 


Log E [GeV] 
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Figure 12. Optical depth (upper panel) and suppression factor (lower panel) 
for absorption of 7 -rays within the jet as a function of the frequency for the 
models of the high states reported in Fig. \m for 6v = 15° (blue) and 
6v = 20° (red). The vertical yellow line shows the upper energy limit of 
the LAT spectrum. 
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Figure 13. The green and the red lines show the emission from the spine- 
layer model of PKS 0521—36 for the cases 6^ = 15° and Ov = 6 ° as 
observed at 0^ = 4°. Red symbols show the high state, green symbols the 
low state. LAT upper limits are not shown. For comparison, the blue data 
points describe the SED of BL Lac itself. 
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6 CONCLUSIONS 

We investigated the properties of PKS 0521—36, an AGN with un¬ 
certain classification. We report results on multiwavelength obser¬ 
vations of this source obtained from radio through 7 rays by SMA, 
UMRAO, VLA, Swift and Fermi-LAI, mostly during 2008 August- 
2012 August. In addition archival XMM-Newton, Chandra, and 
VLBA observations were analyzed. 

The historical classification of this source as a BL Lac object 
is in contrast with the broad emission lines observed in optical and 
UV, with EW larger than 5 A (rest-frame). The presence of broad 
emission lines may suggest a classification as an FSRQ, but the 
radio spectrum below 1.4 GHz is not flat (Or 0.7). The radio 
spectral shape and the optical emission lines indicate a possible 
classification as an intermediate object between BLRG and SSRQ. 
On pc-scales PKS 0521—36 shows a one-sided jet with a brightness 
distribution along the axis decreasing down to the noise level and 
then rising in an extended, stationary feature at nearly 30 pc. This 
is similar to what was observed in M87 and 3C 120, two 7 -ray 
emitting misaligned AGN. The core dominance estimated for PKS 
0521—36 is similar to those of the 7 -ray emitting SSRQ 3C 380 
and 3C 207 and BLRG 3C 111, suggesting an intermediate viewing 
angle with respect to the observer. 

The fit of the XMM-Newton EPIC pn spectrum confirmed the 
presence of a contributi on of the thermal component identified by 
iBirkinshaw et alJ ( l2002h in the Chandra image as a diffuse halo 
around the AGN. No significant Fe Ka emission was detected in 
the XMM-Newton spectrum. This seems to disagree with the possi¬ 
ble classification of PKS 0521—36 as a BLRG, although a smaller 
viewing angle with respect to the other BLRG may lead to a sig¬ 
nificant increase of the jet component that overwhelms the Fe Ka 
line in this object (similarl y to the case of the narrow -line Seyfert 1 
galaxy PMN J0948-l-0022: lD’Ammando et alj2015l) . 

The contemporaneous flux increase observed from optical to 

7 rays during 2010 June suggests the identification of the 7 -ray 
source with PKS 0521—36. The 230 GHz light curve showed an in¬ 
crease of activity peaking on 2011 January. A delay of ~6 months 
between the emission at mm and 7 rays is unlikely, suggesting that 
different parts of the jet are responsible for those activities. How¬ 
ever, we cannot rule out that a mm peak contemporaneous to the 
7 -ray one was missed due to the lack of observations at 230 GHz in 
2010 May-June. The average 7 -ray photon index of the source (T-, 
= 2.45 ± 0.02) is similar to the average value observed for FSRQ as 
well as for the FRII radio galaxies (in particular the SSRQ). More¬ 
over, the average apparent isotropic luminosity of the source (L-, 
= 5.2 X 10^^ erg s“^) is in agreement with the values observed for 
BL Lacs and radio galaxies. However, the lower luminosity with 
respect to SSRQ may be due to the close proximity of the source. 
Thus, taking into account its low redshift, the 7 -ray properties of 
PKS 0521—36 are compatible with those observed in SSRQ and 
BLRG. In this context the strong flaring activity detected from this 
source by Fermf-LAT starting from 2010 June, with daily peak flux 
of l- 2 x 10 “® ph cm“^ s~^, is very intriguing. 

We discuss the 7 -ray emission of this source in the frame¬ 
work of the ‘spine-layer’ scenario, with a fast spine surrounded 
by a slower layer, like in the case of the radio galaxy NGC 1275 
jTavecchio & Ghisellinil2014h . We compare the SED of the flaring 
state (2010 June) with that of a low activity state (2010 February- 
March). We present three alternative models, corresponding to 
three different choices of the viewing angles 9^ = 6°, 15°, and 20 °. 
For the case with 9^ = 6 ° and 15° we obtain a good fit. According 
to the unification scheme for radio galaxies and blazars, reporting 


the SED obtained with 9„ = 15° as recorded by an observer at 9^ 
= 4°, this SED would correspond to that of the ‘aligned’ coun¬ 
terpart of the radio galaxy. On the contrary, the resulting SED is 
strongly unbalanced toward the synchrotron luminosity (by a factor 
of « 100), at odds with the typical blazar SED. This suggests that 
a viewing angle between 6 ° and 15° is preferred. In the spine-layer 
scenario proposed here the rapid variability observed during 7 -ray 
flares favours a relatively small angle. However, we cannot rule out 
the possibility that PKS 0521—36 is an AGN seen at 9v >^15° and 
its ‘aligned’ counterpart is a blazar belonging to a population of 
synchrotron-dominated objects. 

This intermediate viewing angle is in agreement with the core 
dominance and the other radio properties observed, and thus with 
the classification of PKS 0521—36 as an SSRQ at low redshift 
or a BLRG. The strong 7 -ray flaring activity observed by Fermi- 
LAT from this source indicates that SSRQ and BLRG may have 
relativistic jets as powerful as blazars. We noted that an in- 
crease of the 7 -ray acti vity was observed for the SSRQ 3C 380 
jTorresi & Grandiir2013h . although at a flux level lower than that 
reached by PKS 0521-36. 

The multifrequency observations and SED modeling pre¬ 
sented here give important information about the characteristics of 
the AGN PKS 0521—36 and its classification. Further radio-to- 7 - 
ray observations will be fundamental to investigate in even more 
detail the nature and the physical processes occurring in this source. 
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